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Significant reduction in target strength and radiation signature can be achieved by surrounding an 
object with multiple concentric layers comprised of three acoustic fluids. The idea is to make a finely 
layered shell with the thickness of each layer defined by a unique transformation rule. The shell has 
the effect of steering incident acoustic energy around the structure, and conversely, reducing the 
radiation strength. The overall effectiveness and the precise form of the layering depends upon the 
densities and compressibilities of the three fluids. Nearly optimal results are obtained if one fluid has 
density equal to the background fluid, while the other two densities are much greater and much less 
than the background values. Optimal choices for the compressibilities are also found. Simulations in 
2D and 3D illustrate effectiveness of the three fluid shell. The limited range of acoustic metafluids 
that are possible using only two fluid constituents is also discussed. 



o 



C/5 



T3 

C 

O 

o 



> 

O 
O 

q 

o 
o 



X 
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I. INTRODUCTION 

The idea behind transformation acoustics is that a co- 
ordinate transformation makes it possible to have one 
region of an acoustic fluid mimic another region. Fluids 
that have this property have been called acoustic metaflu- 
ids. In transformation optics the transformation uniquely 
defines the material properties, but this is not the case 
in acoustics, and there is an added degree of freedom in 
the makeup of the acoustic metafluid. The range of pos- 
sible acoustic metafluids has been derived^, and includes 
fluids with anisotropic inertia and pentamode materials. 

Interest in transformation acoustics has been moti- 
vated by the possibility of acoustic cloaking. The first 
electromagnetic wave cloaking device^ uses transforma- 
tion of coordinates in the governing wave equation to 
steer energy around the cloaked object. It was sub- 
sequently demonstrated that the same methods should 
work for the acoustic wave equation". The acoustic 
cloak corresponds to the limiting case of a point trans- 
formed into a finite region, and it has unavoidable physi- 
cal singularities associated with the extreme nature of the 
transformation. Different types of singularities are ob- 
tained depending on whether the transformed metafluid 
is purely inertial with anisotropic density and a single 
bulk modulus, or in the other limit, purely pentamodal 
with isotropic inertia. The distinction is important for 
cloaking, for which it is known that use of only flu- 
ids with anisotropic inertia (inertial cloaks) requires in- 
finite mass, and is therefore not a realistic path towards 
acoustic cloaking^. Despite this limitation, it is possible 
to achieve almost perfect, or near-cloaking, using layers 
of anisotropic fluids that approximate the transformed 
medium, without the singularity. For instance, Torrent 
and Sanchez-Dehesapi partition the shell into many small 
but equally thin layers where the local properties are de- 
fined by two normal fluids, with density and bulk mod- 



uli {pj,Kj}, 
tities p r = 



j = 1,2, such that the averaged quan- 
lOi + P2), P± = [UPi 1 + P2 -1 )]" 1 and 



a ' Electronic address: norris@rutgers.edu 



K = [\{K^ + i-CT 1 )] -1 yield the anisotropic metafluid 
properties {p r (r), p±{r), K(r)} proposed by Cummer and 
Schurigi. In order to achieve this equivalence it is nec- 
essary to make {pj,Kj}, j = 1,2, functions of r, with 
the result a large number of distinct fluids is necessary: 
100 and 400 for the two numerical examples reported by 
Torrent and Sanchcz-Dehcsa^. 



The purpose of this paper is to demonstrate that signif- 
icant reduction in target strength can be achieved using 
layers comprised of only three acoustic fluids. The idea 
is to make a finely layered shell that surrounds the struc- 
ture, with each layer being one of the three fluids, but 
instead of prescribing the relative thickness of each layer 
we allow it to be a function of r. The transformation 
formulas then imply unique values for the relative con- 
centrations as functions of r, in both two (cylinder) and 
three (sphere) dimensions. 



The outline of the paper is as follows. The homoge- 
nized layered shell and the transformation metamaterial 
are introduced separately in Section |TT] in the context of 
an N— fluid material. The remainder of the paper con- 
centrates on the 3-fluid (N = 3) configuration. General 
results for both cylindrical and spherical shells are de- 
rived in Section [ill! including the unique transformation 
formulae. Dependence of the cylindrical transformation 
metamaterial on the constituent properties of the 3-fluids 
is explored in Section IIVI The explicit nature of the 
transformation formulae for 2D suggest optimal choices 
for the fluid densities and compressibilities. These find- 
ings are confirmed in Section [V] where examples of cylin- 
drical and spherical 3-fluid nrctamaterials are presented. 
Numerical simulations showing their effectiveness in re- 
ducing scattering strength in 2 and 3 dimensions are also 
presented in Section IVl 
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II. PRELIMINARIES 

We consider radially symmetric configurations, cylin- 
drical in 2D and spherical in 3D. A fluid annulus or 
shell occupies < ro < r < r out , and is surrounded 
by a uniform acoustic medium with density and sound 
speed Pout, c ou t, in r > r out . The shell is assumed to be 
made of a finite number, N, of distinct fluids arranged 
in a well defined stratification that results in an effective 
material with smoothly varying properties in the radial 
direction. We are particularly interested in finding the 
smallest number N for which it is possible that the strat- 
ification has the properties of an acoustic metafluid. An 
acoustic metafluid is defined here as a material with de- 
sirable effective properties that cannot easily be obtained 
with a single, physical fluid. This definition obviously in- 
cludes materials obtained by a coordinate transformation 
of a larger region of uniform acoustic fluid with properties 
equal to those of the exterior fluid in r > r out . 

For simplicity, but with no lack in generality, we set 
Tout = 1, Cout = 1 and p ou t = 1, which is equivalent to 
choosing units for length, time and mass, respectively. 
For the remainder of the paper all quantities are non- 
dimensional. 

We first consider the homogenized shell composed of a 
layering of N distinct fluids defined by their mass den- 
sities, pi, . . . , pn, and the compressibilities C\, . . . , Cn- 
The compressibility is C; = where Ki is the bulk 

modulus, and the wave speeds are Cj = \J ' Ki/ pi, and the 
impedances are Zi = y/Kipl, i = 1, . . . , N. We define, for 
later use, Si = piCi, or alternatively, Si = c~ 2 , so that 
we may identify \/Si as acoustic slowness in fluid i. 

The layering yields an effective fluid with compressibil- 
ity C* and anisotropic inertia defined by radial density 
p r , and circumferential density p±. The parameters of 
the effective fluid arc defined by homogenization of the 
stratified medium asX 



(1) 



where (•) is the local average over the volume fractions 
of the TV-fluids, 





/ (p) 
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It is assumed that 4>i = 4>i(r), so that the averages (fTJ) 
define parameters p r {r), p±(r), and C»(r). This type 
of inhomogeneous or localized homogenization may be 
achieved by allowing the layering to be sufficiently fine, 
and will be illustrated by numerical examples later. 

The transformation from the current (physical) do- 
main to the mimicked one makes the shell appear acous- 
tically as if it is a larger shell of fluid with uniform prop- 
erties equal to the exterior fluid. The key is a transfor- 
mation function, r — > R = R(r), such that the range of R 
exceeds its domain, i.e., the inverse mapping R — > r phys- 
ically contracts space. To be specific, the outer boundary 
is mapped to itself, r = R = 1, and the inner boundary 



r = r is mapped to R = R , with < i?o = R(ro) < r o- 
The perfect acoustic cloak is defined by Rq = 0. The 
transformed material has properties p r x, p±r, and C*r, 
with values uniquely defined by the transformation in 
d— dimensions as^ 



/ Ptt\ 



Rl 



'(r/R) d - ls 
(r/R) 3 ~ d 
Xr/R) 1 -' 



d = 2 or 3, 



(3) 



and where R' = dR/ dr. 

The connection between the homogenized material ([T]) 
and the acoustically transformed material ([3]) is now 
made explicit by requiring p r T = p r , P\_t = P±_ and 
C*t = C* (and we drop the subscript T). Our objective 
is to find families of transformation functions R = i?(r), 
4>i = 0i( r ) f° r which this equivalence can be achieved. 
It depends, of course, on the choices of material proper- 
ties {pi, Ci}, i = 1, . . . , N, and not all combinations will 
work. Among the requirements are that the transforma- 
tion function is one-to-one, and that the volume fractions 
are all between zero and unity. We therefore require that 
<j> £ $at where </> is the N— dimensional vector of volume 
fractions, and $n the N—l dimensional surface on which 
it must lie, 



= {0 4 >O, 5> 4 = 1, i = l,...,N}. 



(4) 



III. THE THREE FLUID MATERIAL 

A. Algebraic formulation 

The first two relations in (TTJ) and the identity @ 2 may 
be written in matrix form for N = 3, 



f 



f 



f 



Pi P2 P3 
.Pi 1 P2 1 P3 1 





(5) 



This can be solved to give the 3-vector of volume frac- 
tions in terms of p r and ■ Substitution into the third 
relation in (fT]) yields an expression for C* in terms of p r 
and pj 1 . Thus, 



^ = /o+Pr/i+P ± 1 / 2J 

C* = a + j3ip r + /? 2 pl\ 

where the 3-vectors in (I6al) are 



(6a) 
(6b) 
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with D = P1P2P3/KP1 - P2XP2 - P3XP3 - Pi)], and the 
scalars a, f3\ and $2 in (]6b[) are 



a=C T f , Px=d F f l , 



c T f 2 



(8) 



with C T = (Ci,C 2 ,C 3 ). 
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B. The transformation function 

1. Differential equations 

Eliminating C*, p r and pj 1 from (|6bl) using the identi- 
ties ([3]) yields a differential equation for the transforma- 
tion function, 



d_R 

dr 



a 



ft-^-ft) 



2D, 
3D, 



(9) 



subject to the boundary condition i?(l) = 1. The pa- 
rameters /? and, for later use, A, /i, are defined 



/3 = ft+ft, A = a + ft /x 



.0 
a 



2. 2D solution 



We first consider the 2D equation (JU^. Let x 
X = R 2 , then eq. ([9]^ becomes 



v dx X 
d X a a 



= 1. 



Integrating yields 



i? 2 + (A - l)i? 2 
A~ 



(10) 



r 2 . 



(11) 



(12) 



The 2D transformation function is therefore completely 
defined by the two parameters A and /i, given in explicit 
form in (IA6I). 



3. 3D solution 

The 3D equation © 2 becomes, with the change of vari- 
able s = 



-as 



1 d.R 

#"dT ~ fas' 1 + as 3 + fas 2 - 1 

= E-^- fl(l) = l., (13) 

where the four roots Sj and the coefficients 7^,1 = 
1, 2, 3, 4, are defined by 

4 

ft Y[(s - sj) = fas 4 + as 3 + fas 2 - 1, (14a) 
— a." 

It 



fa Ilisi-sj)' 



(14b) 



Note that £V 7* = 0, s l = -a/ fa, Y,iH s i = ~a/fa, 
Si7i s ? — (a/A) 2 - Integration of (fl~3f yields 



4 / r , \ 7i 

*4 



(15) 



This provides an implicit formula for R and r, in terms 
of the three parameters a, ft and ft. Using the fact that 
1 < s < so, where sq is defined in the next subsection, 
cq. (|15l) gives R as a function of s, from which r = sR is 
obtained. 




FIG. 1. The range of r/> for the 3-fluid. 



C. The inner radii ro and Ro 

It follows from continuity of the solution of the differ- 
ential equation © that the values of the inner radii ro 
and Rq should correspond to a point on the edge of the 
triangular region $3, see Fig. [TJ The actual radial values 
can be determined from eq. ([5]), using p r and p± as de- 
fined in §3§ and keeping the parameter s = to express 
6j of (l6al) in the form 



p^ 1 + p jPk s 3 - d )R' - ( Pj + p k )] 

(Pi ~ Pj)(Pi ~ Pk) 



(16) 



where i ^ j ^ k ^ i. Replacing R' by © and setting 
(|16p to zero implies an algebraic (polynomial) equation 
for s. In principle there are three possible solutions, cor- 
responding to each of ^ = 0, i = 1,2,3. However, in 
practice for a given set of 3- fluids only one is important, 
and we choose the 3-fluid properties so that it is the root 
for 02 = 0. We consider first d = 2. 

In the 2D cylindrical configuration the equation 02 = 
is a quadratic in s with a single positive root greater than 
unity (corresponding to ro > i?o), which combined with 
(|12p implies Rq and ro in explicit form as 



Ro 



'0 



(A-l) 



Pr2 




- P 


1 - 


Pr2 






1 


P^ 




2 



1-/J 



Ro, 



where 



Pj + Pk 
1 + Pjpk 



(17a) 
(17b) 

(18) 
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For the 3D spherical case the equation 4>2 = becomes 
a biquadratic in s. We find 



i=l 



s - Si 



r = sqRq, 



(19) 



where sq is a positive root of 

s 4 [a+ /3i(pi+ p 3 )]+s 2 [ap 1 p 3 + (3 2 (pi + p 3 )}- (pi + p 3 ) = 0. 

(20) 

The transformation requires So > 1, and numerical ex- 
periments (see Section [Vj indicate that a single real root 
greater than unity exists in all the cases considered. 

D. Total mass and average density 

The total mass m of the 3-fluid shell is the integral 
of the local average of the density, (p). Therefore, m 
follows from eq. ([1]) as the volumetric integral of p r {r). 
Substituting from §5§ i and using ^ , the integral can be 
expressed in closed form for the 2D case, and reduced to 
an integral in s = r/R for the 3D case. We find 



a , (A-i) 



4 / s _ s .f.-' (21) 
4^/s 6 II 1 V ds, 3D, [ > 

1 i=l (i-S'J 

from which the average density in the shell, p = 
3m/[Tr(d + 1)(1 - r#)], can be found. For 2D we find, 
after some simplification, 



(22) 



E. Summary 

We have shown that the three fluid shell is uniquely re- 
lated to possible transformation functions in both 2- and 
3-dimcnsions. The connection is still somewhat tentative, 
since we must confirm that the functions arc physically 
realistic. This requires among other things that the vol- 
ume fractions are all positive and between zero and unity, 
i.e. that <f> £ $3 where the equilateral triangle surface $3 
is defined by (QJ. We must also confirm that the inner 
radii are actually given by eqs. (JTTJ) in 2D and ([TT)| in 
3D. Optimally, both of the inner radii should be small, 
since Ro <S 1 means that the mapped region Ro < R < 1 
is almost the entire interior of the cylinder/sphere of ra- 
dius 1, while Ro <C ro implies that the shell ro < r < 1 
in physical space occupies a relatively small proportion 
of the mapped region. In the next Section we consider 
the 2D shell for which these questions can be answered 
in explicit form. 

The results for the 3-fluid shell indicate that there are 
no free parameters for a given set of fluids. This suggests 
that the transformation property cannot be achieved 
with only two fluids. It is shown in Appendix B that the 
2-fluid case is too constrained, although it does display 
some interesting physical properties, even if it cannot 
provide acoustic cloaking. 



IV. THE THREE FLUID MATERIAL IN 2D 

A. Range of material parameters 

The relation p r p± = 1, which holds only in 2D (sec cq. 
©), considerably simplifies the algebra of the problem 
as compared with the 3D case, allowing clearer under- 
standing of parametric dependence. We refer the reader 
to Appendix A for the details and provide only the main 
findings here. 

With no loss in generality, see Appendix A, we assume 



Pi > P2 > ps, with pi > 1, p 3 < 1. 



(23) 



The density with the intermediate value, pi, may be less 
than, equal to, or greater than unity. In order to dis- 
tinguish these two cases without being specific as to the 
particular one, we define p p as the density with value on 
the same side of unity as p2- We assume for the mo- 
ment that P2 ^ 1; the special case of p 2 — 1 is discussed 
separately below. 

The main result is that the physically obtainable ma- 
terial properties can be parameterized in terms of the 
radial density p r , which has a well defined range itself. 
Thus, <j> £ <& 3 for p rp < p r < p r 2, where p r i are defined 
in (fT5| . The lower bound is not achieved in practice 
but is instead set by the value of p r at r = 1, see Ap- 
pendix A. The physically reachable values of the volume 
fractions, compressibility and density p± are therefore 
defined through p r as 



<t> = 4>rp 



P± = Pr 1 > 



Pr Prp 

Pr2 — Prp 

Pr Prp 

Pr2 ~ Prp 
rv 

for 



)(4>r2 ~<l>rp)> 
) (C*2 — C* p ), 



1-/3 



< 0„ < 



Pr S Pr2, 



(24a) 
(24b) 
(24c) 



where the critical values p ri of the radial density are de- 
fined in (TT5|) , and the critical values of the concentrations 
4>, and compressibilities C«, are 



<Pri = 



>p r =p 



c*\ Pr=Pri =c T <i> r 



(25) 



Based on the sensitivity analysis in the Appendix, the 
radial density p r has its greatest range, defined by Ap r = 
Pr2~ Prp, if Pi is large, P2 is close to unity, and p 3 is small. 
Thus, 



Pi > P2 ~ 1 > P3 



1 + P1P3 



(26) 



The optimal strategy seems to have three fluids with 
properties in line with (f26|) . For instance, if (pi, p2,Ps) = 
(10, 1.1, 0.01), then Ap r = 8.002 as compared to 8.182 ac- 
cording to (f2"6")l . The scalings of (|2^) also imply that the 
relative magnitudes of the concentrations of the light and 
heavy fluids are 



03 (r) « pip 3 <f>i(r). 
It is also shown in the Appendix that 

(C*2 — Pr2)(C*p — Prp) < 0, 



(27) 
(28) 



4 




FIG. 2. The range of $ for the 3-fluid in the cylindrical con- 
figuration. The dashed lines show the possible straight line 
paths as a function of p r . In practice, the path begins at some 
point inside the triangular region (r = R = 1) and ends at 
<j> 2 =0 (r = r ,R = Ro). 



which places another constraint on the choice of the three 
fluids involving their compressibilities. We next consider 
these extra degrees of freedom in the context of a special 
case of (|26[) which makes the parameterization simpler. 



B. The case of p 2 = 1 and other limits 

The previous results, in particular the suggested opti- 
mal strategy for choosing the densities of the three fluids, 
suggests that the results will not depend strongly on p 2 
if it is close to unity. It is therefore reasonable to simply 
take p 2 = 1, which leads to other simplifications which 
we now examine. 

The reachable line in $3 has one end at the vertex 
62 = 1, and (pMl) becomes 



(29a) 
(29b) 















4>r2- 
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r 







K Pr2 - 1' 



c 2 + (c* 2 -c 2 )(J^-), 

Pr2 — I 



for the same range of p r as in , and with 



4> 



\i-pDpi 

•2 ~ — 2 2 I ^ 

P^-Pa V(p?-iw 



Pr2 



C*2 




The two parameters in the transformation function ([12 
simplify, using, (|A6|) . to 



A = 5*2, fi = l- 



(Pi - - - P3)5 2 



(l- /3 2 )|Sl + (/9 2_ 1)53 _ (p 2_^ )S . 2! 

where 5; = piCj, i = 1,2,3. Equations (JTTJ) simplify for 



p 2 = 1, using (|A12|) , to give 
S 2 



i?0 



so = 



1 



S§S 2 - 1 



r = sqRq, 



(1 - pftS! + {p\ - 1)S 3 



9 2 

Pi - Pa 



(30a) 
(30b) 



We next examine these exact results for some limiting 
values of the other 3-fluid parameters. 

Both quantities in (|30|) should be small. Based on the 

assumed density scalings (f2"Bl) . it follows that — §■ can be 

r o 

small only if both Si =0(1) and S3 =o(l). Under these 

circumstances, eq. (|28)) . which is now (S 2 — l) (l — > 

0, requires S2 > 1, and ([30|>9 implies in turn that /i ~ 0. 
We therefore have, in addition to (|26[) for the densities, 
that the quantities Si, i — 1,2,3, should satisfy Si = 
0(1), S 2 > 1 and S 3 < 1. 



1. The case p 2 — 1, pi — p 3 1 

Further simplification results from setting = p^ 1 , 
still with pi 3> 1. For instance, the volume fractions of 
phases 1 and 3 are equal, 

h=h = ki-h) = k^— \), (3i) 

Z Z p r2 — 1 



and p r 2 reduces to p r2 = \ (pi + p x 1 ) ■ 



2. Summary 

Based on the analysis above it appears that optimal 
choices for the properties of the three fluids are 

Pi > 92 = 1 > Ps, Si = 0(1), S 2 > 1, S 3 < 1, 

(32) 

implying A = S2, p ~ 0. Under these circumstances, 
(|30p provides the relatively simple approximations for the 
values of the inner radius ro, and its pre-transformed 
value, Rq, 



ro 



1 



(33a) 



« {(! - irK 5 3 + S lP f)} 1/2 « 1. (33b) 

The value of ro can be made to be close to unity by 
further requiring 

S 2 > 1 => r « 1-^-, Ro w (Sa + Sipr 2 ) 172 - (34) 

ZD2 

For this range of parameters the thickness of the physical 
shell, 1 — ro ~ , depends only on the squared slowness 
S2, while the image of the inner radius, Ro, is dependent 
on the other two slownesses, and the density p\. While 
the parameters ro and Ro are insensitive to the densities 
p 2 =0(1) and P3 =o(l), the other quantities, such as C» 
and <j) can depend on these. However, if pz = 1/pi then 
the concentrations of fluids 1 and 3 are everywhere the 
same. 
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TABLE I. The four cases of 3-fiuid material considered. 



Case 


Pi 


pi 


PS 


Si 


s 2 


5*3 


1 


10 


1 


0.2 


1 


10 


0.1 


2 


10 


1 


0.2 


1 


10 


0.01 


3 


100 


1 


0.02 


1 


10 


0.01 


4 


1000 


1 


0.002 


1 


10 


0.01 



V. NUMERICAL RESULTS 



A. Example of three-fluid shells 




The range of possibilities for the 3-fluid metamaterials 
is extensive given that there are 3x2 = 6 independent 
variables at our disposal. However, based on the esti- 
mates in Section IIV1 particularly (|32[) , it seems reason- 
able to take p2 = Si = 1. We further take p3 = 2/p%, 
in keeping with (|32j) . Also, considering f|34[) we choose 
S2 = 10, which leaves two parameters: pi and 5*3. Four 
distinct 3-fluids are considered according to the four sets 
of parameters in Table I with different combinations of 
pi and S3. The transformation functions and the con- 
centrations of the three fluid constituents are illustrated 
in Figs. [3]- [6] The curves R = R(r) illustrate the trans- 
formation, which maps the original region Rq < R < 1 
to the physical domain ro < r < 1 , and the values of the 
inner radii, ro and Ro, are given in Table II. Note that 
R < r, as expected. Also, the concentrations for the 2D 
shells, in Figs. [5^ and H^i, satisfy $3 « 2</>i, in 

accordance with (|27[) since P1P3 = 2. The most impor- 
tant aspect is the relative values of ro and Rq, in that it 
is desirable to have ro close to unity while Ro should be 
close to zero. The value of ro is smallest in Fig. [3] and 
largest in Fig. [6l and it appears to increase with p\ . In 
order to obtain a value of ro close to unity, and in good 
approximation with the estimate p4[b . it is necessary to 
have a large value of pi, see Figs. [5] and [6] Although 
only two values of S3 arc considered here, numerical ex- 
periments indicate that the value of R is more sensitive 
to this parameter, with Ro decreasing as 5*3 is increased. 
It is also found that better results, i.e. smaller Ro, larger 
ro, are obtained when S2 becomes very large. For in- 
stance, ro = 0.989, Ro = 0.031 is obtained in 2D with 

pi = s 2 = 10 3 , s 3 = 10- 3 . 



1 

0.8 
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FIG. 3. (color online) The curves show the concentrations of 
the three fluids and the radius R as functions of the physical 
radial coordinate r for the fluid parameters of Case 1 (see 
Table 1). (a) the 2D cylindrical configuration; (b) the 3D 
spherical shell. 



TABLE II. Results for the four cases of Table I. p is the aver- 
age density in the shell ro < r < 1. uo is the relative value of 
the total scattering cross section at kro = 3 of a rigid cylin- 
der/sphere surrounded by the 3-fiuid shell with 500 layers. A 
value of 100% corresponds to the bare rigid target. 







2D 








3D 








ro 


Ro 


P 


ffo(%) 


ro 


Ro 


P 


<ro(%) 


1 


0.60 


0.20 


3.12 


25.8 


0.66 


0.26 


5.41 


4.55 


2 


0.41 


0.06 


3.13 


2.37 


0.59 


0.19 


5.69 


2.20 


3 


0.88 


0.09 


19.17 


0.69 


0.88 


0.11 


57.7 


.033 


4 


0.94 


0.09 


40.22 


0.69 


0.96 


0.096 


192 


.012 



B. Discrete layering algorithm 

The inhomogeneous nature of the homogenized mate- 
rial is captured by layering the shell on two scales. The 
first scale is a fine layering of L distinct bands defined by 
the regions between ro < ri < r 2 < ■ . ■ < = r out = 1. 
The second scale of layering defines three sub-regions be- 
tween neighboring radii. Let r n \ = r n , and define 

r i, m = r t, m -i - <t>m-i{r n ) A n , m = 2, 3, (35a) 
A„ = r^-r^_ 1 , n=l,2,...,L, (35b) 

where j(d+ 1)A„ is the area or volume between the in- 
ner and outer radii of the band [r n _i, r n ]. The three re- 



gions (r n . 2 ,r n s] , (r„, 3 ,r„, 2 ] and (^-1,1,^,3] have frac- 
tional volumes <pi{r n ), (f>2(r n ) and 4>3(r n ) of the band, 
respectively, and are therefore occupied by the respec- 
tive fluids, see Fig. [7] The choice of the ordered set 
{r„, n = 1,2, . . . , L — 1} is relatively arbitrary as long as 
it is finely spaced for large values of L. For simplicity 
we take A n constant, independent of n, in which case 
A„ = (1 — Tq)/L = A and the radii become 

r t,i = r o + nA , n=l,2,...,L, (36a) 
ri, m = Tn,m-1 ~ 0m-l(»n,l)A, m = 2, 3. (36b) 
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C. Scattering from a three-fluid shell 

We consider plane wave incidence in the uniform ex- 
terior fluid r > 1, with time harmonic dependence e~ lult 
(henceforth omitted). The 3-fluid shell in tq < r < 1 
is defined by the discrete layering algorithm, and is as- 
sumed to surround a rigid object of radius ?*o. The scat- 
tered pressure is expressed 

oo 

p(r,9) = J2 A nMr)K(0), (37) 

71=0 

where 9 is the polar angle with respect to the incident 
direction, {ip n (r), A n (9)} = {Hn\kr), cos n9} in 2D and 

{h i n ) (kr),P n (cos9)} in 3D, and k = to is the nondimen- 
sional wavenumbcr. In the shell region the pressure p 
and radial velocity v are expressed in modal form 

oo 

(p(r,6), v(r,9)) = £ (p n (r), v n (r))A n (9), r < r < 1. 

n=0 

The 2- vector U(r) = {p n {r), r d ~ 1 v n (r)^ satisfies the 
ordinary differential equation (ODE) 

*E=Q(r)U(r) t r <r<l, (39) 
d r 

where 



and % 2 = ?i 2 , ti(ti + 1) in 2D and 3D, respectively. The 
density p(r) and bulk modulus K(r) are piecewise con- 
stant, defined by the 3-fluid material properties at each 
value of r according to the discrete layering algorithm. 



1. Computational scheme 

Three different numerical methods are employed to 
find the scattered pressure (|37|) : (i) by solving for the 
matricant; (ii) using a global matrix; and (hi) by solv- 
ing the matricant of the homogenized radially dependent 
anisotropic fluid. In the first method the matricanl£, or 
propagator matrix, is found by numerical integration of 
the matrix equation dM/dr = QM subject to the ini- 
tial condition M(ro) = I, the 2x2 identity matrix. Then 
using the continuity conditions at r = 1, and the rigid 
boundary conditions at r = 7*0, it is possible to express 
the scattering coefficient A n in terms of M (1). Solution 
(ii) using the global matrix method, e.gi^, is obtained by 
creating a large system of simultaneous equations which 
can be cast as a matrix equation of size 6L. The third 
method (ii) is based on the equations of motion of an 
anisotropic acoustic fluid, e.g£, with radially varying pa- 
rameters p ri p± and C* given by the exact transforma- 
tion formulas ([3]). The equations of motion can be trans- 
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FIG. 8. (color online) Case 1. The magnitude of the scattered 
pressure for an incident wave of unit amplitude for the 2D 
(top) and 3D (bottom) 3-fluid shells. In each case kro — 3 
and L = 500. The inner dark circular region depicts the 
rigid target of radius ro, surrounded by the shell of unit outer 
radius. 



2. Numerical results 



FIG. 7. The discrete layering algorithm to reproduce the local 
homogenization properties of the 3-fluid shell. 



formed into the form (|39|) with Q — > where 



QJr) 



P d-1 



„2d-4/ 



(41) 



and = C^ 1 . Using an ODE solver it is again possible 
to find the scattering coefficient A n . Details of numerical 
schemes (i) and (iii) will be provided in a forthcoming 
paper. 



Figures [8] and [9] show the magnitude of the scattered 
acoustic field for an incident wave of unit amplitude. 
Since the radius of the object being cloaked changes for 
each of the four cases of Table I we take the nondimen- 
sional characteristic value kro = 3 in each scattering sim- 
ulation. This allows us to compare the total scattering 
cross-section between the four cases even though the val- 
ues of tq are different. Fig. [10] shows the response of 
the bare 3D spherical rigid target based upon case 3 in 
which ro = .88. The total scattering cross-section for 
the "cloaked" rigid object was calculated using the coef- 
ficients An, and compared with the cross-section for the 
bare rigid object. In each Table II shows, the 

relative cross-section is diminished, and for cases 3 and 
4 the reduction is significant. Note that the reduction in 
target strength is greater in 3D as compared with 2D, in 
agreement with the general findings o&. The numerical 
methods (i) and (ii) were found to be in agreement with 
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FIG. 9. (color online) Case 3. The same as for Figure [8] 2D 
and 3D simulations are in the upper and lower plots, respec- 
tively. 



one another, and with method (iii) when L is very large. 
For instance, the cross-section found using method (iii) 
is 0.3% larger than that of method (i) for the 2D exam- 
ple in Fig. [5] Finally Fig. [TT] shows the effect of the 
number of layers L on the relative value of the total scat- 
tering cross section for case 3. A curve fit of the power 
function aL b + c shows that for this particular case at 
this particular frequency the cross-section decreases as 
do ~ L~ 2 ' 2 . More layers provide a better approximation 
to the homogenized limit, as expected. For small num- 
bers of layers the layering algorithm used here could be 
improved using various optimization strategies, but we 
do not pursue that here. 

VI. CONCLUSION 

The main finding of this paper is that it is possible to 
achieve cloaking-likc behavior with as few as three dis- 
tinct acoustic fluids. Using transformation acoustics, we 
find that for a given set of three fluids the layered shell 
ro < r < 1 is uniquely determined, with the inner radius 
ro given by eqs. (fTTj) and ([!§]) for 2D and 3D, respec- 
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FIG. 10. (color online) 3D pressure map solution for a rigid 
cylinder; kro = 3, ro = .88. 
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FIG. 11. (color online) The relative value of the total scat- 
tering cross section for case 3 vs. the number of three- 
fluid layers in ro < r < 1. The dashed lines are curve 
fits of the form aL b + c for 5 < L < 100. The pa- 
rameters (a, b, c) and the root mean squared error (RMSE) 
were found to be (3716,-2.221,0.9924), 0.290 for 2D, and 
(6435, -2.258,0.1324), 0.278 for 3D. 



tively. The shell is made of fine layers of the three fluids 
with relative concentrations as a function of r determined 
from eqs. (HJ), (|6a|) and (|9|). Obviously, the overall ef- 
fectiveness and the precise form of the layering depends 
upon the relative densities and compressibilities of the 
three fluids. The best results are obtained if one fluid 
has density equal or close to the background or host fluid 
density, while the other two densities are much greater 
and much less than the background value. Numerical 
simulations of the scattering from specific layering real- 
izations confirm the theoretical predictions and show the 
effect of the finite number of layers. Many questions re- 
main as to the the optimal choice of fluids in general, and 
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what can be achieved using existing fluids specifically. 
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APPENDIX A: PROPERTIES OF THE 2D 3-FLUID 
MATERIAL 

1. Density and compressibility 

We begin with the density implications. Equation ([fjal 
reduces, using p r p± = 1, to give 



Pi(Pj + Pk)(Pr ~ Pri) 
Pri(pi ~ Pj)(Pi ~ Pk) ' 



i + 3 ^ k ^ i, 



(Al) 



where the critical values of p r are given by (|18p . Based 
upon the identities (|A1|) . we note that 



0, 



u 3\p r =c 



Pj Pri 



Pi 



p\ 



(A2) 



where i =^ j =^ k ^= i. The points defined by (|A2|) are 
the intersections of the line (|A1|) with the planes • (j> = 
0. In order to have some <j> £ $3 at least one of the 
intersections must lie on the boundary of $3. Consider 
p r i of (|18[) . then <pj and <pk must both be positive, which 
occurs if and only if one of (pj,Pk) is larger than, and 
the other is less than, unity. This gives an important 
necessary condition: At least one of the three densities 
is larger than unity, and conversely, at least one must 
be less than unity. This condition must hold in addition 
to the obvious requirement that the three densities are 
distinct, since otherwise the system ([5]) is not solvable. 

Introduce the density values p p , p m , {p 7^ m} <G {1, 3}, 
such that 

(pi - l)(Pp - 1) > 0, (p 2 - l){ P m - 1) < 0, (A3) 

with 2^p^m^ 2. 

We note some other properties of the critical values of 
the densities: 



Pri Prj — PriPrj 



1 



( Pj - Pi ){i - pD 
(pi + pk)(pj + PkY 

{Prj -!)(!- Prk) 
1 ) 
Prj 1 Prk 



(A4a) 
(A4b) 



where i =/= j =/= k =/= i. These imply, respectively, that 

Pri > Prp > Prm, and p r 2 > 1, Prp > 1 and p rm < 

1. Combining these with the previous inequalities, we 
surmise the ordering p r 2 > Prp > 1 > prm- Thus, for 
instance, if p2 > 1, then the possible range of p r is p r \ < 
p r < p r 2- If pi < 1 then it is p r z < p r < Pri- 

Any value of p r in the range p rp < p r < Pri therefore 
yields a triple of concentration values satisfying ^ £ $3. 
At the upper (lower) value, p r = p r i (= p r p)i the concen- 
tration <f> lies on the boundary of the triangle with <f>2 = 



(4>p = 0). But these limiting values are not necessarily 
achieved. Thus, at r = R = 1 the differential equality 
© implies that p r = a/(l - (3), see eq. §M\i . This is 
the practical lower bound on the range of p r . Equations 
(|24|) and ([25]) 1 then follow. 

By analogy with equation for the volume fractions, 
the effective compressibility of (|6b[) can be expressed in 
the form (|24b[) . Alternatively, eq. (|6b[) implies C* = 
a + Pp r , and therefore we deduce that a and (3 may be 
expressed 



PrlC* p PrpCt, 1 C* 2 C* p 

a = , \i = . 



Pri — prp 



Pri Prp 



(A5) 



These lead in turn to explicit expressions for the two pa- 
rameters that define the transformation function, (|12[) . 



A _ 5i (l-Pa)(l-ft») , s ^ (l-pa)(l-Pi) 



{Pi- Pl){Pl- Pz) (p2 -Ps){p1-Pl) 

(l- Pl )(l-p 2 ) 



(A6) 



3 (P3- Pl)iP3- Pi)' 
PrlSl{pl - P§) + PrlHpl - Pi) + PrzS^pl - p\) 



sM-pl) + Si{pi-p{) + Sz{pi-pi) 

where Si = piCi, i = 1,2, 3. 
a. Sensitivity 

The reachable range of p r is, from (J9j) , p rp < p r < 
p r p + Ap r where 

Apr = pri- p rp = ^ , „ „ W1 , „ „ ■ (A7) 



(1 + pip m )(l + PpPrn) 



Hence, 

dApr (1 - p 2 J 



dpi (1 + pip m ) 2 

OApr (1 - P 2 J 



dp p (1 -I- PpPm) 2 ' 

(9Ap r _ (p p - p 2 )(l + P 2 m )(p2 +Pp + ZpiPpPrn) 



(A8) 



dpr, 



(1 + P2/3m) 2 (l + PpPm) 



If p = 1 these are, respectively, < 0, > 0, < 0. Con- 
versely, if p = 1 they arc > 0, < 0, > 0. Hence, whether 
p = 1 or p = 3 it is clear that Ap r is greatest if pi is 
large, pi is close to unity, and P3 is small. 



2. Transformation function 

a. Necessary conditions for the three-fluid parameters 

Since i?(l) = 1 at the outer radius r = 1, we have 

i?'(l) = Pr (l)=a(l) = T ^, (A9) 

that is, 

Pr (l) = C,(l) - PnC.r-prrCrt _ (A1Q) 
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But we require that p rp < p r (l) < Pr2, or, since p r2 
Prp > 0, 



C*p Prp 



p r (l) - p rp _ 

Pr2 Prp Pr2 Prp + C*p C* 2 
Pr2 — Pr(l) _ Pr2 — C*2 

Pr2 ~ Prp Pr2 — Prp + C» p — C» 2 

Hence, ([25]) must hold, or, explicitly 



>0, 
> 0. 



(All) 



£3— t~ 

Pl - P3 Pi - P3 

!-Pm , c P^ 1 



1 x 



Pi - Pm 



Pi - Pm 



1) < 0. 



(A12) 



b. The case p 2 = 1 

In this case the p/m distinction is unnecessary since 

Prl = Pr3 = 1 for P2 = 1, =^ G $3 for 1 < /9 r < /5 r2 . 

This implies that the reachable quantities reduce to 
(tig) . 

APPENDIX B: THE TWO FLUID MATERIAL 

1. General theory 

The 2-fiuid version of eq. (|5|) is 



(Bl) 




This implies that the concentrations are 



P2 



Pi - P2 



Pi - Pr 
Pi ~ 92 ' 



(B2) 



and the densities p r , p±_ arc related by the compatibility 
condition for (IBll). 



(B3) 



Pr + PlP2P± = Pi + P2- 



The effective compressibility, which follows from (|B2[) 
and the third relation in (fTJ), satisfies 

(pi - P2)C* + (C 2 - d)p r = Pl C 2 - paCi. (B4) 



Equation (|B2|) provides relations for the volume fractions 
in terms of the radial inertia p r . One can also interpret 
cqs. (|B3[) and (|B4[) as defining pj_ and C*, respectively, 
in terms of p r . Therefore, all parameters in the two- fluid 
material can be defined by a single quantity, in this case 

Pr- 

However, in order to relate the two-fluid material to 
a transformation it is necessary that there exists a func- 
tion R which satisfies the three differential identities . 
Substitution of these into eqs. (|B3[) and (|B4[) gives a 
pair of equations which can be considered as algebraic 
equations in two unknowns: R' and R/r. Solutions for 



both of these quantities can be found in terms of the two- 
fluid properties pi, p2, C\, C%, but the solutions are not 
of practical interest. The reason is that the constant val- 
ues of R' and R/r that arc found, say R' = a, R/r = b, 
must be equal, leading to trivial cases. The main conclu- 
sion from the study of the N = 2 case is that the 2-fluid 
material is overly restrictive. 



2. A special case of a uniform 2-fluid material 

While it is not possible for the 2-fluid material to re- 
produce a transformation material, it is possible to make 
some interesting uniform fluids with anisotropic inertia. 
The idea is to seek constant values of p r , p± and C* which 
also match to the exterior fluid in r > 1. This requires 
that R = 1 at r = 1. Enforcement of © then requires 
the three parameters in the left vector be equal to R'. 
Substituting into eqs. (|B3[) yields 



= Pj 1 = C*= p r i 



The volume fractions follow from (IB2I) as 



PlPrZ 



\-p\ 



Pi 



p\ 



P2pr3 



P\ ~ Pi 



(B5) 



(B6) 



which arc both positive if and only if (1 — pi)(l — p 2 ) < 0. 
The one remaining condition, for the compressibility, im- 
plies using (|B4j) and ()B5j) that the two compressibilities 
must be related such that 



C lPl (l - pi) + C 2 p 2 (p\ -\) = p\- pi 



(B7) 



The anisotropic fluid (|B5[) is defined by the parameter 
Pr3 = Prs{pi, P2), and is composed of volume fractions 
4>i = 4>i{pi,P2) of fluid i = 1,2. Denote any pair sat- 
isfying the relation (|B7|) as d = Ci(pi,p 2 ), i = 1,2. 
It is interesting to note that these functions are in- 
variant under the interchange {pi, p 2 , <pi, 2 , C\, C 2 } — > 



{P2 



<Pl 1 ! 



</>2,</>l,C 2 ,Cl}. 



a. Examples 

If, for instance, C\p\ = C 2 p 2 then (|B7[) implies that 
C\p\ = C2P2 = !■ Both fluids have the same wave 
speed as the background fluid. They differ only in their 
impedances, which in this case are = pi = , 
i = 1,2. 

Conversely, if C\/ p\ = C 2 /p 2 then (|B7[) implies that 
Ci/pi = C 2 1 p 2 = 1 . The two fluids have the same acous- 
tic impedance as the background fluid, and differ only in 
their wave speeds, which are Ci = p~[ = , i = 1,2. 



3. A two and a half fluid material 

As a case intermediate between the strictly 2-fluid and 
3-fluid cases, consider the 2D case, for which (r/R)R' = 
p r , see J3Ji. It follows from (|B3[) 9. i.e. p r = , that 
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Pr = Pr3, a constant. Taking into account the boundary 
condition R(l) = 1, the unique mapping is 

R(r) = r Pr3 . (B8) 

Equation (|B4|1 combined with ©3 then implies 

(1 - pDS, + {p\ - l)S 2 = [p\ piy^z-V. (B9) 

This cannot be satisfied if the two fluids have properties 
independent of r. However, if we still require that the 
densities are fixed, but the compressibilities could vary 
with r, then (|B9|) suggests that a mapping can be realized 
if one or both Si, £2 are such that the equality holds for 
some range of r. It is well known that adding a small 
concentration of bubbles to a liquid results in an increase 
in the compressibility without significant change in the 
effective density. Hence, it might be possible, in principle 
if not in practice, to add a third fluid whose only role is 
to enhance compressibility. In this sense it is half of a 
fluid, since its inertial properties are not used. 
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